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INTRODUCTION 


During  the  period  covered  by  this  report  four  distinct 
types  of  related  projects  were  in  progress  in  this 
laboratory , involving  the  interaction  of  electromag- 
netic radiation  with  plasmas.  They  were: 

(a)  Electron  temperature  variations  in  a current  modu- 
lated plasma  column. 

(b)  Afterglow  electron  density  decay  mechanisms  and 
energy  losses  in  high  energy  argon  and  helium  plasmas, 
including  measurements  of  ion  conversion  rates. 

(c)  Computational  study  of  highly  ionized  plasmas  in- 
volving a one-dimensional  code  for  the  simulation  of 
hydrogen  and  helium  afterglows. 

(d)  Generation  and  detection  of  far  infrared  (sub- 
millimetre radiation  using  a HCN  gas  laser  and  semi- 
conductor detectors. 

Since  most  of  this  work  has  already  been  or  is  in  the 
course  of  publication  in  journals  or  conference  pro- 
ceedings, this  report  will  summarise  the  research  by 
means  of  abstracts  of  the  papers,  augmented  by  more 
detailed  descriptions  where  necessary. 
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THE  VARIATION  OF  ELECTRON  TEMPERATURE  IN  PLASMA  COLUMNS 


The  method  involved  the  design  of  a sensitive  11.5  GHz 
microwave  radiometer  and  is  fully  described  in  the  follow- 
ing papers: 

1)  An  electronically  switched  microwave  radiometer 
6. J. Burk ley,  C.P. Downing  and  M.C. Sexton 
J.Phys.D  (Appl.Phys),  6,  668,  1973 

Abstract: 

An  11.5  GHz  microwave  radiometer  is  described  in  which 
periodic  noise  variations  up  to  1 kHz  were  resolved  by  a 
sampling  process  using  PIN  diode  switching. A sensitivity 
of  0.04  dB  (100 was  recorded  in  checking  the  system 
with  the  noise  power  emanating  from  a current-modulated 
plasma  column. 

2)  Electron  temperature  variations  in  deeply  modulated 
plasma  columns 

C. J. Burk ley , and  M.C. Sexton 

Bulletin  of  Am.Phys.Soc.  18,807,1973 

(Presented  at  25th  Annual  Gaseous  Electronics  Conf., 

London, Ontario, 17-20  October  1972) 

Abstract : 

An  11.5  GHz  microwave  radiometer  was  used  to  monitor 
electron  temperature  variations  in  argon  and  krypton 
plasma  columns  by  means  of  a sampling  process  involving 
PIN  diode  switching.  Temperature  variations  as  low  as 
O.Ci  -IE  (100°K)  were  recorded  in  0.5A  columns  current 
modulated  up  to  85%  with  sinusoidal  and  triangular  wave- 
forms. 

3)  Electron  temperature  measurements  in  a current 
modulated  plasma  column 

C.J.Burkley,  D.M. Barry  and  M.C. Sexton 

Proc.llth  Int.Conf.on  ion  phenomena  in  gases, Prague, 

September  1973,  p.42. 

Abstract : 

The  measurement  of  electron  temperature  in  the  positive 
column  of  a gas  discharge  was  carried  out  by  monitoring 
the  microwave  noise  radiation  with  a sensitive  11.5  GHz 
radiometer  using  Schottky  barrier  mixer  diodes.  PIN 
diode  switching  was  introduced  to  enable  the  noise  rad- 
iation from  the  discharge  to  be  compared  with  a cali- 
brated noise  source  under  current  modulated  (sine  and 
triangular)  conditions  in  the  positive  column. 

In  argon  and  krypton  plasmas  and  with  mean  currents  of 
approximately  0.5A,the  modulation  was  shown  to  produce 
a "hysteresis"  effect  in  the  current/temperature  relation. 
This  was  ascribed  to  differing  electron  radial  distribu- 
tion for  increasing  and  decreasing  modulation  currents. 
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A feature  of  this  work  (which  has  now  terminated  - 
see  Interim  Scientific  Report  No.2,p.l)  is  that  although 
the  radiometer  performance  was  checked  with  a modulated 
plasma  only, it  can  be  applied  to  any  system  in  which 
rapid  power  fluctuations  are  likely  to  occur,  eg., 
fading  on  microwave  radio  links. 


AFTERGLOW  DECAY  MECHANISMS 


(a)  Electron  density  measurements  and  afterglow 
thermal  economy: 

The  design  and  construction  of  a Hydrogen  Cyanide 
(HCN)  laser  which  operates  at  890  GHz  (337yM)  and  a 
submillimetre  interferometer  with  an  8 millimetre 
microwave  bridge  which  together  measure  electron 
densities  in  the  range  7 x 1014  cm'3  to  8 x 101Acm”3 
in  the  afterglow  of  pulsed  helium  and  argon  discharges 
was  described.  Electron  temperature  measurements  were 
made  on  both  gases  with  a monochromator -photomultiplier 
system. 

Recombination  coefficients  were  extracted  from  the 
electron  density  decay  curves  and  these  indicated  the 
dominance  of  collisional  radiative  recombination  in 
both  gases.  In  helium  the  electron  temperature 
dropped  from  70O0oK  to  less  than  3000°K  in  0.6  msec, 
while  the  recombination  coefficient  varied  from  almost 
10“9  cm3  sec"1  to  7 x-lo"'11  cm3  sec*"1  over  the  electron 
density  range  3 x lOi*  cm”3  to  5 x 1019  cm”3.  A com- 
parison of  these  results  with  previous  experimental 
work  on  helium  showed  good  agreement.  The  effect  of 
including  an  elevated  and  time-varying  neutral  gas 
temperature  in  the  calculations  involving  ion-atom 
energy  exchange  brought  theoretical  predictions  and  the 
present  experimental  results  into  excellent  agreement. 

The  measured  values  of  electron  temperature  in  argon 
were  much  higher  than  in  helium, dropping  from  over 
16,000  °K  to  8,000  °K  in  0.35  msec.  The  experimentally 
determined  values  of  recombination  coefficient  varied 
from  1.6  x 10"9  cnr  sec”*  to  1.54  x 10"11  cm3  sec”1 
over  the  electron  density  range  8.25  x 1011  cm”3  to 
3.6  x 1014  cm“3.  These  results  compared  very  favour- 
ably with  previous  experimentally  determined  values. 

In  the  calculations, the  effects  of  optical  transparency 
to  plasma  radiation  were  also  included. 

Papers ; 

4)  Afterglow  decay  processes  in  highly  ionized 
helium'  and  argon  plasmas 
J.H. Mount joy  and  M.C. Sexton 
Bull . Am. Phys . Soc.  19,159,1974 
(Presented  at  26th  Gaseous  Electronics  Conf., 
Madison, Wis.  October  1973) 

Abstract: 

A HCN  laser  interferometer2was  used  to  monitor  electron 
densities  from  1015  to  lO1^  cm""J  in  He  and  Ar  after- 
glows in  which  electron  and  neutral  gas  temperatures 
were  also  monitored.  Binary  electron-atom  (collisional) 
recombination  was  clearly  identified  as  the  principal 
decay  mechanism, with  coefficients  within  the  range 
10"  to  lo”  cm3  sec”1.  Excellent  agreement  with  theory 
was  obtained  by  including  elevated  and  time-varying  gas 
temperatures  into  the  calculations. 
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(b)  Ion  conversion  in  high  density  afterglows 

A free  space  microwave  interferometer  was  used  to  moni- 
tor the  decaying  electron  density  afterglows  of  argon  and 
helium  plasmas  pulsed  by  varying  d.c. voltages . The  formation 
of  the  molecular  ion  Ar^  was  found  to  depend  on  the  input 
power  - a result  not  previously  reported  - which  should  have 
significance  for  developments  in  high  power  rare  gas  lasers. 

In  the  present  work  the  experimental  parameters  were 
as  follows: 

Microwave  probing  frequency:  35  GHz  (cut-off  density 

cm" ^ 

Maximum  discharge  conditions  20kV,  5kA  for  20ys 
(variable) : 

Discharge  tube  diameter:  10.2  cm  (internal) 

argon  at  0.3  - 1 torr 
helium  at  1 - 6 torr 

The  basic  reaction  for  atomic  to  molecular  ion  conversion 

during  an  afterglow  is  , 

X+X+X^Xj+X 

where  the  second  atom  X,  i.e.  the  third  body, carries  off  the 
excess  energy  and  accordingly  stabilises  the  reaction. 
Clearly  the  reaction  is  a function  of  pressure  and  since 
electron  recombination  to  atomic  ions  (X  j is  considerably 
less  than  that  to  molecular  ions  (X*)  a rapid  decrease  in 
the  afterglow  electron  decay  time  occurs  if  the  reaction 
occurs  during  increase  of  neutral  gas  pressure. 

Paper 

5)  Ion  conversion  rates  in  the  afterglow  of  high  power 
argon  miasmas 

P.J. Murphy  and  M.C. Sexton 
Bull.Am.Phys.Soc.  20,256,1975 

(Presented  at  27th  Gaseous  Electronics  Conf ., Houston, 
Texas,  October  1974) 

The  dependence  of  the  reaction  rate  constant  k for  the  pro- 
cess Ar+  + 2Ar  Art  + Ar  was  examined  with  a 35  GHz  micro- 
wave  interferometer  over  the  electron  density  range  lo!3  - 
loll  cm"  3 in  argon  afterglows.  The  value  of  k varied  from 
2.9  to  1.5  x 10" 31  cm6  sec"l  as  the  input  power  increased 
from  15  to  350  mJ  cm"3.  Moreover, the  gas  pressure  at  which 
the  onset  of  the  conversion  process  occurred  Increased 
during  the  same  range  of  input  power. 
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Similar  results  in  helium  gave  k decreasing  from 
3.5  - 1.4  x 10“32  cm6  sec-1  as  the  electron  temperature 
varied  from  350  °K  - 1300  °K. 

Both  the  afterglow  and  ion  conversion  work  de- 
scribed in  this  section  have  now  been  concluded. 


COMPUTATIONAL  STUDY  OF  HIGHLY- IONIZED  PLASMA  AFTERGLOWS 

Throughout  the  experimental  afterglow  work  described 
above,  computer  simulation  studies  (Glowcode)  were  ini- 
tiated in  collaboration  with  the  Culham  Laboratory 
(Euratom/UKAEA  Fusion  Association)  and  Oxford  Polytechnic. 
The  GLOWCODE  program  simulates  the  behaviour  of  an  after- 
glow plasma  in  one-dimensional  cylindrical  geometry  and 
its  purpose  is  to  study  and  clarify  the  various  competing 
processes  taking  place  in  an  afterglow  in  which  inertial 
effects  may  be  neglected.  The  numerical  model  predicts 
the  radial  and  temporal  distributions  of  a three  fluid 
plasma  comprised  of  positive  ions, neutral  atoms  and  elec- 
trons. Thus  the  afterglow  is  described  by  the  dependent 
variables  ne,ni,na,  Te»Ti»Ta  which  represent  the  electron, 
ion  and  atom  densities  and  the  electron, ion  and  atom 


temperatures  respectively. 

For  computational  purposes  the  plasma  consists  of  a 
nesting  set  of  cylindrical  shells, each  shell  character- 
ised by  its  mean  radius  and  its  thickness .The  calculation 
proceeds  by  calculating  the  energy  changes  due  to  the 
various  collisional  and  radiative  processes  in  each  shell, 
and  the  transport  between  neighbouring  shells.  At  each 
stage  of  the  calculation  pressure  balance  is  restored  by 
adjustment  of  the  cylindrical  shells  until  the  total 
particle  pressure  gradient  is  zero  everywhere.  Since 
the  plasma  is  subject  to  a continual  radial  motion  it  is 
convenient  to  express  the  various  physical  equations  in 
Lagrangian  form. 

The  code  is  structured  within  the  OLYMPUS  package 


introduced  by  Christiansen  and  Roberts  (Com.  Phys.Comm. 
7,245,1974).  Using  this  system  it  is  possible  for  the 
user  to  make  ad  hoc  changes  to  the  program  using  the 
EXPERT  facility  so  that, for  example, different  rate  co- 
efficients could  be  prescribed  for  a particular  physical 
process.  Alternatively, a particular  physical  process  may 
be  suppressed, which  is  useful  when  one  is  studying  the 
relative  importance  of  the  various  competing  processes. 
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Paper;  GLOWCODE:  a one  dimensional  code  for  the 
simulation  of  plasma  afterglows 

J.W.Long,  A. A. Newton  and  M.C. Sexton 

Com. Phys. Conun.  9,1976  (accepted  for  publication) 

This  is  a comprehensive  paper, the  main  features  of  which 
are  abstracted  above.  The  computational  details  are  as 
follows. 

Title  of  Program : GLOWCODE 

Program  obtainable  from;  CPC  Program  Library, Queen's 
“ university  of  Belfast, N. Ireland 

Computer:  ICL  System  4:  Installation » Culham  Laboratory 

Operating  System:  Multijob 

Programming  language  used : STANDARD  FORTRAN 

No. of  bits  in  a word:  32 

High  Speed  Store  required:  24000  words 

Overlay  Structure:  None 

No, of  Magnetic  Tapes  required:  None 

Other  peripherals  used:  Line  printer 

No. of  cards  in  combined  program  and  test  deck:  2487 

Card  punching  code:  EBCDIC 

CPC  Library  Subprograms  used:  Catalogue  No;  Title 

ABUF  OLYMPUS 

Ref. in  CPC: 

7(1574)245 

Keywords : Afterglow, ionisation , recombination , diffusion, 
implicit , Lagrangian . 

Following  the  setting  up  of  the  program  theore- 
tical studies  were  carried  out  on  hydrogen  and  helium 
afterglows. 

7)  Computational  study  of  highly  ionized  plasma  after- 
glows: 

J.W.Long,  A. A. Newton  and  M.C. Sexton 
1st  Europhysics  Study  Conf.on  Atomic  and  Molecular 
Physics  of  Ionized  Gases,  Versailles , April  1973, 
Paper  F6. 

Results  are  presented  for  the  hydrogen  afterglow  with  the 

following  parameters: 

Electron  density  : 10^**  - 10*^  cm  ^ 

Afterglow  period  : 0 - 100  ysec 

Gas  pressure  : 0.2  - 20  torr 


8. 


The  significance  of  recombination  and  atom  thermal  con- 
duction is  confirmed, with  radiation  losses  contributing 
greatly  in  the  early  afterglow.  Electron  loss  by  dif- 
fusion and  wall  recombination  is  small  and  the  observed 
electron  density  gradient  is  due  to  diffusion  of  heat 
rather  than  of  particles.  Above  1 torr  the  incidence  of 
molecular  ion  formation  - especially  significant  in 
hydrogen  - has  also  been  included. 

8)  Numerical  investigation  of  highly  ionized  helium 
plasmas 

J.W.Long,  A. A. Newton  and  M.C. Sexton 

2nd  Europhysics  Study  Conf.on  Atomic  and  Molecular 

Physics  of  Ionized  Gases,Innsbrflck,Sept.l974,p.22. 

The  decay  processes  in  a highly  ionized  helium  plasma 
have  also  been  studied  numerically  in  one-dimensional 
cylindrical  geometry  for  starting  temperatures  of  the 
order  10,000  °K  with  gas  pressures  20-500  mtorr  and  the 
same  density  range  as  for  hydrogen.  Experimental  work 
had  established  that  the  decay  is  dominated  by  three- 
body  collisional-radiative  recombination  In  which  the 
ionization  energy  is  transferred  to  the  electrons;  sub- 
sequently this  energy  is  lost  by  atom  thermal  conduction 
to  the  containing  wall.  The  situation  in  helium  is 
complicated  by  the  possibility  of  recombination  to  one 
of  the  helium  metastable  states, in  particular  the  2^S 
level.  The  helium  plasma  is  assumed  to  contain  electrons, 
neutral  atoms, singly  and  doubly  ionized  ions. Radiative 
transfer  is  not  considered  but  radiation  loss  processes 
are  included.  As  for  hydrogen, the  main  experimental 
decay  parameters  were  confirmed  especially  the  effect  of 
the  2^S  metastables  in  maintaining  elevated  electron 
temperatures  during  the  afterglow. 

9. 


■ 

GENERATION  AND  DETECTION  OF  SUBMILLIMETRE  RADIATION 

As  referred  to  on  p.4  of  this  report, the  main  plasma  diag-  . 

nostic  has  been  the  hydrogen  cyanide  (HCN)  gas  laser 

interferometer  operating  at  337  microns.  This  enabled 

15  -3 

electron  densities  up  to  approximately  5 x 10  cm  to 
be  readily  monitored  in  the  highly  ionized  plasmas. 

In  addition,  as  a result  of  direct  contact  with  the  French 
Commissariat  a l’Energie  Atomique,  a special  HCN  laser 
system  was  designed  and  constructed  for  use  as  a diagnostic 
for  thermonuclear  machine  plasmas. 

Since  1975,  attention  has  been  concentrated  in  Cork  on  the 
laser/detector  system  itself.  The  laser  was  pulsed  in  the 
region  of  10  kV,  10  - 100A  to  ascertain  input  powers  suit- 
able for  optimum  lasing  and, in  addition, the  mean  radial 
electron  density  (10*1  - I010cm3)  during  lasing  was 
measured  with  an  8 mm  microwave  interferometer.  (It  will 
be  noted  that  the  pulse  network  and  microwave  equipment 
previously  developed  for  the  plasma  afterglow  work  have 
now  been  adapted  to  monitor  the  laser  plasma  itself) . 

Mathematical  models  of  the  laser  amplifying  medium,  i.e. , 
the  plasma, were  formulated  to  calculate  parameters  such 
as  preferred  oscillation  modes,  mode  compression  near  the 
critical  density  and  onset  of  instability  with  particular 
reference  to  both  radial  and  longitudinal  electron  density 
distributions.  A comprehensive  analysis  is  in  course  of 
preparation, but  some  preliminary  results  have  already  been 
presented  in  Europe  and  in  the  USA. 

9)  Gain  and  Visualization  of  the  modes  of  a Thermally 
Stabilised  HCN  Laser 
J.P.Lesieur,  M.C. Sexton  and  D.Veron 
J.Phys.D, Applied  Physics,  5,1212,1972 

Abstract:  Optimization  of  a 337  ym  HCN  laser  involving 

long-term  stability , choice  of  beam  splitter  (polyethylene) 
and  maximum  power  has  been  carried  out.  The  gain  deter- 
mination leads  to  a value  of  4%  per  metre.  In  addition, 
a liquid  crystal  camera  has  been  designed  to  demonstrate 
visually  the  characteristics  of  several  inodes  of  the 
cavity. 
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Papers ; 


10)  Mode  characteristics  of  a pulsed  HCN  laser 
X.Hamonic  and  M.C. Sexton 

Proc.l2th  Int.Conf.on  Ionization  Phenomena  in  Gases, 
Eindhoven, August  1975, p. 221. 

Abstract;  A study  of  the  line  emission  from  a pulsed  HCN 
laser  revealed  a strong  dependence  on  the  repetition 
frequency  which  was  attributed  to  the  continuing  presence 
of  excited  HCN  molecules  up  to  20  milliseconds  after  the 
current  pulse.  An  "inversion  density"  time-constant  of 
approximately  2 milliseconds  was  deduced. ..Microwave 
measurements  of  the  electron  density  OfclO1^  cm-3)  during 
the  lasing  action  enhanced  the  explanation  of  the  ob- 
served spectrum  and  demonstrated  the  existence  of  mode 
compression  near  the  critical  density. 

11)  Influence  of  the  electron  density  on  the  modes  of 
a pulsed  submillimetre  laser 

A.Hamonic  and  M.C. Sexton 
Bull.Am.Phys.Soc.  21,142,1976 

Abstract;  The  line  emission  at  337pm  from  a pulsed  HCN 
laser  (%1  torr,5kV,25A  for  25psec)  was  monitored  with 
pyroelectric  detectors.  Simultaneous  electron  density 
measurements  during  the  lasing  action  were  made  with  an 
8 mm  interferometer.  Both  the  asymmetry  of  the  line 
shape  and  cyclic  variations  in  output  power  with  pulse 
frequency  strongly  indicated  the  presence  of  excited  HCN 
molecules  up  to  20  ms  into  the  afterglow.  Moreover, 
calculations  taking  radial  electron  distributions  into 
account  supported  quantitatively  the  presence  of  mode  com- 
pression near  the  critical  density. 


The  availability  of  suitably  fast  detectors  of  sub- 
millimetre radiation  is  of  paramount  importance.  Pyro- 
electric detectors  have  been  successfully  used  in  the 
present  work.  However, during  1975/76  a photoconductive 
detector/FET  preamplifier  unit  at  cryogenic  temperatures 
was  successfully  developed  in  association  with  the  Culham 
Laboratory.  Both  gallium  doped  and  boron  doped  germanium 
crystals  have  recently  been  tested  at  66  microns  wave- 
length. Design  criteria  include  (a)  bandwidth  up  to  50 
MHz  approximately  (b)  noise  reduction  and  minimization  of 
external  capacitances.  This  work  is  now  nearing  comple- 
tion and  will  be  fully  reported  in  due  course. 
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